Entropy-Driven Phase Transitions of Entanglement 
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We study the behavior of bipartite entanglement at fixed von Neumann entropy. We look at the 
distribution of the entanglement spectrum, that is the eigenvalues of the reduced density matrix 
of a quantum system in a pure state. We report the presence of two continuous phase transitions, 
characterized by different entanglement spectra, which are deformations of classical eigenvalue dis- 
tributions. 



PACS numbers: 03.67.Mn, 02.50.Sk, 68.35.Rh 

Entanglement is an important resource in quantum in- 
formation and computation [1]. For bipartite systems, it 
can be quantified in terms of several physically equiva- 
lent measures, such as purity and von Neumann entropy 
[2]. To optimize the use of this resource, towards quan- 
tum applications, it is important to understand which 
states have large entanglement and how these states can 
be produced in practice. A proper understanding of ran- 
dom pure states is crucial in this context. Random pure 
states are known to be characterized by a large entan- 
glement and a number of interesting results have been 
obtained during the last few years. 

Lubkin [3] understood that the average purity of a bi- 
partite system is almost maximal if the pure state of the 
total system is randomly sampled. The analysis was ex- 
tended to higher moments by Giraud [4] , and to the aver- 
age von Neumann entropy by Page [5]. All these results 
are a consequence of a fundamental phenomenon, namely 
the concentration of the entanglement spectrum, that is 
the eigenvalues of the reduced density matrix [6]. The 
typical entanglement spectrum at fixed purity was deter- 
mined in [7] , where the presence of some phase transitions 
was unveiled. This result was extended to different Renyi 
entropies by Nadal et al. [8]. 

In this Letter we further extend these findings to the 
von Neumann entropy. This is the final step towards a 
thorough understanding of the typical bipartite entangle- 
ment of pure states. We shall see that this step, besides 
having a fundamental interest, also discloses results that 
are somewhat unexpected. We shall find two phase tran- 
sitions. Remarkably, one of them, related to the "evap- 
oration" of the largest eigenvalue, is softer than in the 
case of purity (and all other Renyi entropies), in that it 
becomes continuous. 

Our calculation hinges upon the saddle point equations 
for a partition function and makes use of a Coulomb gas 
method [9]. It is valid when both subsystems are large. 



A byproduct of our results is the probability distribu- 
tion of the von Neumann entropy of random states, or in 
other words, the relative volumes of the manifolds with 
constant entanglement (isoentropic manifolds). 

We consider a bipartite system in the Hilbert space 
W = Ha <8> Ha> described by the pure state \ip). The 
reduced density matrix of subsystem A is the Hermitian, 
positive, and unit-trace matrix 

A = tr^|V>><# (1) 

The bipartite entanglement between A and A is quanti- 
fied by the von Neumann entropy of qa 

N 

■S'vn(A) = -tr(gAlngA) = - ^ A fc lnA fc , (2) 

fe=i 

where N = dim Ha, A = (Ai,...,Ajv) G Ajv-i are 
the Schmidt coefficients (i.e. the eigenvalues of qa), and 
Ajv_i is the simplex of eigenvalues (A& > 0, J2k = !)• 
We are interested in balanced bipartitions: TV = 
dimV.A = dim Notice that < SVn < In A, where 
the minimum and maximum values are obtained, re- 
spectively, for separable and maximally entangled vector 
states. 

We will focus our attention on the typical properties of 
the eigenvalues A of qa- For random pure vector states 
sampled uniformly on the unit sphere (ip\ip) = 1, the 
eigenvalues of the reduced density matrix are distributed 
according to the joint probability density function [10-12] 

Pjv(A) = C N [] (Xj - A fe ) 2 , (3) 

l<j<k<N 

Cn being a normalization factor. Starting from (3), Page 
found that the average value of the von Neumann entropy 
is almost maximal, namely, for large N [5, 13] 

S vN = lnA-i (4) 
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Indeed, for large N, the distribution p^ concentrates 
around a typical A, that maximizes pn [6] , and the typical 
spectral distribution of qa for large N is known to follow 
a Marcenko-Pastur law [7, 14] with support [0,4/A]. 

A natural and more general question is how the en- 
tanglement spectrum is typically distributed in a sys- 
tem with a certain amount of bipartite entanglement. 
In other words, one is interested in the typical distribu- 
tion of the Schmidt coefficients on isoentropy manifolds, 
conditioned at a given value of the entropy SVn ■ 

This is a constrained maximization problem: given a 
value u € [0,lniV] find A such that 

PnW — maxpTv with 5 v n(A) = In N — u. (5) 

By introducing two Lagrange multipliers £ and j3, that 
constrain the eigenvalue normalization and the deviation 
u of the von Neumann entropy from its maximum In N, 
respectively, the problem is translated into the (uncon- 
strained) minimization of the potential 



v{U,p) = -J^Eha, - A *i +?(E A * - 1 

j<k ^ k 

+ p(jTx k \nNX k -u^, (6) 



with respect to A, £, and (3. This is the energy of a 
gas of charges (eigenvalues) distributed in the interval 
[0, 1] with a 2D (logarithmic) Coulomb repulsion, sub- 
ject to two external electric fields proportional to £ and 
j3. The logarithmic form of the interaction is a direct con- 
sequence of the product form (3) of the joint probability 
density. 

It is worth noting that this problem can be equivalcntly 
framed in the statistical mechanics given by the partition 
function [7] 



J N 



-L 



-0N 2 h(\) 



Pjv(A) d N X, 



(7) 



with an "energy density" h(X) = InN — 5 v n(A) and an 
inverse "temperature" /?. In the thermodynamic limit 
N — > oo, one looks at the maximum of the integrand, 
that is at the minimum of the potential (6). Large values 
of f3 yield highly entangled states, while (3 = yields 
random states. 

The saddle-point equations, dV/d/3 = dV/d£ = 
dV/d\ k = 0, read 



EVn(A^) 



(8) 



/3(lnA/A fc + l) + A E ' A ^_ + ^ , (9) 

for 1 < k < N, where the primed sum is restricted to j ^ 
k. When all eigenvalues Xj are of order 0(1/N), we can 



introduce the empirical distribution of the eigenvalues 

^^E^-^)' ( 10 ) 

that in the limit of large N can be approximated by a 
continuous probability density function. By making use 
of (10), Eqs. (8)-(9) read 

J A In A cr(A) dA = u, J Act(A) dA = 1, (11) 

/3(lnA+l) + 2^^1dA'+e = 0, (12) 

with A = A^Afe, and / denoting the Cauchy principal 
value. Equation (12) can be solved using a theorem by 
Tricomi [15]. The solution lies in a compact interval [a, b] 
with < a < b, and takes the form 



a(A) 



f r(X'] 

L 7T(A'- 



r(A')_ (b-X')(X'-a) 
A)V (6-A)(A-a) 



dA', (13) 



where r(A) = -(£ + /? + /? In A) /2tt. 

From the second equation in (11) and the conditions 
of regularity at a and b, that are equivalent to 

a(a) = 0, (7(b) = 0, (14) 

the edges of the distribution are readily found to be 
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Moreover, the spectral density reads 



1 



.(A) = ^p^-AKA-,),^ - M),(16) 



1 . (15) 



7r(6-a) 2 
where r\ = (b + a)/(b — a) , and 

rj + ^7] 2 - 1 I 1 ln(y + rj) 



g(x,v) 



2tt 



V 1 - y 2 (y - x ) 



dy (17) 



is a universal function, with x € [— 1, 1] and T] > 1, that 
deforms Wigner's semicircle law (obtained for purity [7]). 
See Fig. 1. 

Notice that the lower end a of the eigenvalue distri- 
bution is positive for (3 > (3 C = 3/2 and vanishes when 
(3 = j3 c . This is a critical value, at which we encounter 
the first phase transition. Indeed, for (3 < (3 C , one gets 



o = 0, b= -(V2/3 + 1-1), 



(18) 



so that the lower end stays still at a = and the eigen- 
value distribution is no longer regular at a = 0. One 
gets 



a(A) 



irb\ X 



(19) 
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FIG. 1. (Color online) Deformation function g in (17) for a 
few values of rj. The function g is defined in (20). 
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FIG. 2. (Color online) /3 vs it, from (22). 

with a deformation function 



g(x) = 2(x + l)g(x,l)-l. 



(20) 



See Fig. 1. In particular, notice that for /3 = 0, b = 4, 
the deformation in the spectral density (19) vanishes, and 
one recovers the classical Marcenko-Pastur law 



cmp(A) 



1 

2^ 



A 



A 



(21) 



The von Neumann entropy at the typical entanglement 
spectrum (5) has a distance u from its maximal value 
that is given, in the limit TV — > oo, by the first equation 
in (11). By plugging (16) and (19) into it one gets 



In 1 



u((3) = 



In 



1 

7 + 1 



1 

2/3 



1 



0</3< 2 



(22) 



where 7 = y/1 + 2/3. The inverse of this function j3 = 
0(u) is plotted in Fig. 2 and enables us to express every- 
thing in terms of the amount of bipartite entanglement as 
measured by the von Neumann entropy SVn = In N — u. 
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FIG. 3. (Color online) Entanglement spectra <r(A) for various 
values of von Neumann entropy S v n = In N — u: (a) < u < 
u c ~ 0.26, (b) u c < u < 1/2, and (c) 1/2 < u < InN. 



In particular, we notice the maximal value u(0) — 1/2, 
which is the average value (4), and the critical value 



3\ 2 2 
".• = - I = In - + - ~ 0.26, 



(23) 



at which the entanglement spectrum changes its phys- 
iognomy, through a continuous phase transition. The 
entanglement spectrum er(A) is displayed for a few values 
of u in Figs. 3(a) and (b). 

One can extend the analysis to the case u > 1/2, 
towards separable vectors with A ~ (1,0,..., 0). In 
the statistical-mechanics model this would correspond to 
negative "temperatures" (3 < 0. By setting Ai = (J, = 
O(l), for large N, the saddle-point equations (8)-(9) re- 
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duce to 



/zlniV 



-V'- 



N 2 ^ Xj 

j>2 J 



(fc > 2), 



(24) 
(25) 



and £ = /3hiN = j3. By means of the empirical distribu- 
tion 



*(A) = 

they become 
Act(A) dA = 1 



V l-u 



7V- 



J">2 



(26) 



A' -A 



dA' - /3(1 - /i) = 0. (27) 



They are equal to (11)-(12) with /3 = and £ = - 
/i). Therefore, besides the eigenvalue fi, the spectrum of 
the reduced density matrix is made of a sea of eigenvalues 
whose distribution cr(A) = <tmp(A) is given by (21). See 
Fig. 3(c). 

Finally, by evaluating the density function at the typ- 
ical spectrum A one gets 



Pat (A) oc e 



(28) 



where s is the entropy density of the statistical-mechanics 
model (the "entropy of the entropy"): 



N 2 ^ 



\n\N\j - NX k \. 



(29) 



In the limit N — > oo one finds 



s = < 



In 



-In 



1 



1 

4' 



1 P 3 

- +7 -2"2' 



ln(l -n)- 



1 

2' 



< u < u r 



Un < U < 



<u< IniV, 



(30) 



where 7 = yT + 2/3, /3 = /3(m) is the inverse function of 
(22), and [i = u/\nN, according to (24). The logarithm 
of the probability (28) (i.e. the volume) of the isoentropic 
manifolds is plotted in Fig. 4 for TV = 50. Observe that s 
is unbounded from below and the isoentropic manifolds 
shrink to a vanishing volume both at u — (maximally 
entangled states) and u = IniV (separable states). The 
probability that a random state be maximally entangled 
is therefore exponentially suppressed as N — > 00. 

The presence of discontinuities in some derivatives of 
the volume detects the two phase transitions. At u = u c 
there is a continuous phase transition, associated to the 
vanishing of the lower edge of the spectrum, and signaled 
by a discontinuity in the fourth derivative of s, as shown 
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FIG. 4. (Color online) Logarithm of the volume of the isoen- 
tropic manifolds a = N~ 2 InpN vs u — In N~ S v n, for N = 50. 
See (30). The discontinuity of the derivative at u = 1/2 is 
0(l/\nN). 
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FIG. 5. (Color online) The third derivative of s with respect 
to u. The fourth derivative jumps at u = u c . 



in Fig. 5. A second phase transition occurs at u = 1/2, 
due to the split off of the largest eigenvalue 0(1) from 
the others 0(l/N). It is worth noticing that this phase 
transition of the von Neumann entropy is softer than the 
analogous one for purity [7] and other Renyi entropies [8] , 
which are first order. Interestingly, when the Renyi expo- 
nent a becomes smaller than 1, this phase transition dis- 
appears. The von Neumann entropy signals therefore a 
crossover between a violent (first-order) phase transition 
towards separable states, associated to the evaporation 
of the largest eigenvalue for a > 1, and the absence of a 
phase transition when a < 1. This and other interesting 
issues will be discussed in a follow-up article. 
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